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We consider the Bak-Sneppen model near zero dimension where the avalanche exponent r is 
close to 1 and the exponents jj, and a are close to 0. We demonstrate that t — l = fj, — a — 
exp{— /i^^ — 7 + . . .} in this limit, where 7 is the Euler's constant. The avalanche hierarchy equation 
is rewritten in a form that makes possible to find the relation between the critical exponents a and 
H with high accuracy. We obtain new, more precise, values of the critical exponents for the ID and 
2D Bak-Sneppen model and for the ID anisotropic Bak-Sneppen model. 
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Perhaps the most simply formulated model showing 
avalanche behavior is the Bak-Sneppen model |]-|: 
"What could be simpler than replacing some random 
numbers with some other random numbers?" . Never- 
theless, the exact solution of the Bak-Sneppen model is 
unknown even in one dimension. The value of the most 
fundamental quantity, i.e. of the upper critical dimension 
is also still under discussion ^,0. 

The formulation of the model is short indeed. A ran- 
dom number from some distribution 'P(/) is placed 
at an each site i of a lattice. One replaces simultane- 
ously the smallest number fmin of them and the random 
numbers at its nearest neighbour sites by new random 
numbers from the distribution V{f), and afterwards the 
process is repeated. 

Avalanches in the Bak-Sneppen model are defined in 
the following way. The /-avalanche is defined as the se- 
quence of the steps at which fmin remains smaller than 
the given parameter /. (One may find a more detailed 
definition in [^.) 

A very significant step to the understanding of the na- 
ture of the avalanches in this model was made in [ p| by 
S. Maslov who introduced the so called avalanche hierar- 
chy equation for the distribution P{s, f) of /-avalanche 
sizes s (i.e. of temporal durations). From this exact 
equation, one may get the additional relation between 
the critical exponents of the model. Unfortunately, the 
exact solution of the equation is known only for the mean 
field situation. Two first terms of the expansion from the 
mean field solution, i.e. from the upper critical dimen- 
sion, were calculated in iQ. The precision of the results 
obtained by direct numerical integration of the avalanche 
hierarchy equation in its original form [ pd| is only compa- 
rable with the precision of the Monte Carlo simulations 

There is another way to get analytical results. It seems 
natural to start from the lower critical dimension which 
equals zero for the Bak-Sneppen model to make some- 
thing similar to the well known 2 + e expansion. Here, we 
derive from the avalanche hierarchy equation some conve- 
nient relations which enable us find the singular relation 
between the critical exponents near zero dimension and 



get values of the exponents at integer dimensions. Tra- 
ditionally, one relates the exponents r and /i (see 
the definition of these exponents below) . The total curve 
t(/z) with the particular points for the integer dimen- 
sions and the areas of applicability of the approaches of 
|]l2| and ours is depicted in Fig. 1. 

For the distribution P{f) = J > @, the 
avalanche hierarchy equation is of the form pd| 
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= Y,t^P{t,J)P{s-tJ)-s^P{sJ). (1) 
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Here, in the scaling region, gives the average num- 
ber of the distinct sites updating during an avalanche of 
the size s, where — d/Df, d is the dimension of the 
lattice, Df is an avalanche fractal dimension The 
physical meaning of the equation describing the hierar- 
chical nature of avalanches in the Bak-Sneppen model 
is the following. The distribution P(s, /) changes while 
/ grows because of two reasons. First, two consecutive 
avalanches of size t and s — t contribute to the avalanche 
of size s (the second avalanche starts from one of the sites 
changed during the first avalanche that gives the factor 
t^^ in the sum). Second, some avalanches of size s merge 
into a larger avalanche. 

For the problem under consideration, the exponent /i, 
which equals at the lower critical dimension and equals 
1 at the upper critical dimension (simulation in joj and 
|l0| , gave different values for it, due = 8 and 4, corre- 
spondingly), is the given parameter. All other exponents 
are related with /i by Eq. (|l|). Below the threshold /c 
(i.e. in the symmetric phase) the solution of Eq. (^ has 
the following scaling form: 



F(s,/) = s--F(s-(/e-/)) 



(2) 



Eq. (|l|) resembles nonlinear differential equations with 
a peaking regime For such equations, it is possible 
to find both exponents included in Eq. (|^). 

In , it was proposed to search for the Laplace trans- 
form of the distribution P(s, /): 
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p(a,/)=^P(s,/)e-"^ 



(3) 



Inserting Eq. into Eq. one obtains the usual 
relation between the critical exponents 



Then Eq. (|l|) gives 
1 dp{a, /) 



of ^ 



fa 



(4) 



where /da^ denotes the fractional partial derivative 
(/i is certainly noninteger) and the last expression is its 
integral representation. r( ) is the gamma- function. The 
scaling relation for the solution of Eq. (0) is 



(5) 
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II — a 



(6) 



and the following integral-differential equation for the 
scaling function h(x): 



T{1- fi)x 
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h(x) 
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Here h'{x) = dh{x)/dx. 

In ||l2| , Eq. (|^) was used to obtain the expansion from 
the mean field solution but it seems to be incon- 

venient. Let us show that one may transfer it to a purely 
integral form. The following lines demonstrate how that 
integration may be done. 
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Applying dx to the first and the last lines of Eq. (H) and then integrating by parts (one may chose h{x = 0) = 1 
II) we get 



r(l - /i) hih{x) 



dz 



[i-z^-Y 



xh{xz) 



fi — a 



duh{uz) 



and finally we obtain the equation for the scaling function h{x) in the most convenient form: 



h[x) = cxp 



1 



dy 



r(i - /i) 7o [1 - [y/xy/'^Y 
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(9) 



(10) 



(if one does not demand /i(0) = 1, h{x) in the left parts of Eqs. (^ and ( [lO|) is h{x)/h{0)). 

Asymptotic form of h{x) for large x follows from the expansion of Eq. (H) in small a. Below the threshold, h{x) 
has to be 



h{x) ^ x'^i'-"^l"-^l''{cQ + cxx-^l'' + cix-"^!" + ...), 



(11) 



This particular asymptotic behavior fixes the solution of Eq. (10) and the value of a for any given /i. Substituting 
Eq. (|ll|) into Eq. (0), Eq. (1), or Eq. (|l|), one gets the sume rule 



dxh{x) = il^^ r(l - fi). 



(12) 



Note that if h{x, /i, a) is a solution of Eq. (hO[) then ch{cx, fi, a) is also a solution for any constant c. 
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Eqs. ( p^ ) and ( [T^ ) are the set of equations that lead to the scahng function h{x, /i) and cr(/i). Instead of Eq. (|T^), 
one may use equally the condition on the value of the exponent of the asymptote 



[xh' (x) / h{x)]{x co) = 



(13) 



Hence, the problem is reduced to the eigen value problem for the nonlinear equation [T^ . 

Let us study the solution of the system for small /i. The expansion of the solution of Eq. (|ic| ) in x looks as 

1 fj, 



lnh{x) = 1 - - B{a,l- fi) 



r(i-/i) 
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B{3a, 1 - n) 



r(i-M) 



(14) 



where B{ , ) is the beta-function. We shall see that 
the quantity {fi ~ a)/a is the smallest parameter of the 
problem near the lower critical dimension. If one tends 
formally n to 0, the last line of Eq. (nj) tends to 



rr ^— ' n n ' 



and afterwards to 



r(i-M) 



i^h{x) = -^y^( , " , 

^ ' a ^ nn\ VT 1 - u / 

n— 1 ^ ^ ' / / 



(15) 



(16) 



Thus, for small /i, the solution h{x) behaves as the follow- 
ing. For low enough x, the solution very slowly decreases 
from the value ft.(0) — 1 and in some crossover region, 
X l//i, it comes to the asymptotic power tail, Eq. (11). 
We have to stress that three last limit equalities may be 
justified only for small x and that the omitted terms in 
Eq. (|l^ do contribute to the power-law tail. Neverthe- 
less, one may try to estimate cr(/x) for small /i inserting 
Eqs. (15) and ( p^ ) into the sum rule (12). The solution 
of the first of the equations obtained in such a way is 
/i — (T = exp{— //^^ — 27 + 0(/i)}, where 7 — 0.5772 ... is 
the Eulcr's constant, and the solution of the second one 
is of the form 



^-a — T~l — exp{— ^ ^ - 7 + 0(/i)} . 



(17) 



Thus, the dependence is non analytical but even the sec- 
ond term of the expansion can not be defined by such 
estimation. 

In fact we failed to obtain the value of the constant 
analytically. Nevertheless, Eqs. (|l^) and ( p^ ) are very 
convenient for numerics, since iterations of Eq. ( [lO| ) con- 
verge. (One may start, for instance, from the functions 
(H) or {^.) We checked the validity of the relation (0) 



for small fi. The value of the constant in Eq. (jlj) ob- 
tained in such a way is 0.5771(5), i. e. that is the Eulcr's 
constant indeed. 

Solving Eq. ( [lO| ) with the constraint, Eq. ( |T^ ) or 
Eq. (^3|), and the initial condition, one may easily get 
a{i-j,) and r(/i) for any given fi. The values of the ex- 
ponent fi are known from simulation at integer dimen- 
sions with much higher precision than the values of r be- 
cause of the better available statistics [||,^ . Therefore 
we can improve essentially the precision of the known 
value of T. For the ID Bak-Sneppen model, we get 
T = 1.0637(5) + 0.4(/i - 0.4114), where fi = 0.4114(2) 
is the value obtained from the Monte-Carlo simulation 
p3{ . For the 2D Bak-Sneppen model, we obtain r = 
1.229(1) + 0.77(/i - 0.685), where /i = 0.685(5) is the 
value obtained in (The last relations may be used 

to obtain better values of r when the more precise values 
of the exponent fi will be available.) Now the precision 
of T coincides with that one of fi. Note that these val- 
ues are below the values of r previously obtained from 
the simulation, r(lD) = 1.073(3) and t{2D) = 1.245(10) 
JlSf but arc in accordance with the less precise values 
found in p2) by direct numerical solution of the avalance 
hierarchy equation, Eq. (Q). The value of the exponent 
T of the 3D Bak-Sneppen model may be obtained from 
the expansion from the upper critical dimension Jl^ . In 
Fig. 1, we show the curve T(/i) together with the points 
for the integer dimensions and the low-/z asymptote, (|r^) , 
and the expansion from the upper critical dimension p2| . 

Of course, the relation ( p7| ) is valid only for /i <C 1. 
Nevertheless, let us compare the value of t at /i = 0.4114 
obtained from Eq. (p^), r = 1.0494, with the calculated 
above t{1D). One may see that these values are in qual- 
itative agreement. 

In the case of the ID anisotropic Bak-Sneppen model 
(i. e. for the update of the extremal site and only one its. 
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for instance, right neighbour) the exponents a and ^ are 
coupled by the following additional relation: a + fj, = 1 
p6t . Hence, one can find all the exponents of the prob- 
lem. Wc obtained from Eq. ( [lo| ) the following value, 
/i = 0.5779(5). In two different values for /i were 
obtained, fi — 0.58 from Eq. (|^) and /i = 0.588 found in 
another way. The Monte-Carlo simulation made in | |T^ ] 
and [|l^ gave = 0.60(1) and fi = 0.59(3) correspond- 
ingly. Therefore, we had to check our result. For that 
we solved numerically Eq. (^ with the initial condition 
h{0) = 1 and the constraint, Eq. (|l|) or Eq. (|l|). The 
result is /z = 0.5778(5). Thus, the value fi — 0.578 seems 
to be more reliable but the question is still open. 

In summary, we have demonstrated that the simple 
transformation of the avalanche hierarchy equation made 
it convenient for analysis and numerics. We have ob- 
tained the non trivial singular relation, r— l = /i — (t = 
exp{— — 7 -|- . . .} with the Euler's constant 7, be- 
tween the scaling exponents of the Bak-Sneppen model 
near zero dimension. Using the known from simulation 
values of the exponent /i, we have found, in fact, all other 
exponents of the Bak-Sneppen model in ID and 2D with 
the same high precision. We have got also the exponents 
of the anisotropic ID Bak-Sneppen model. 

Nevertheless, one should note that the main problem 
of obtaining of the last independent critical exponent of 
the Bak-Sncppcn model in a regular way remains open. 
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FIG. 1. The exponent r vs. fi calculated from Eq. ( |lC| ) 
(see also [p^J^). The value of the exponent fi depends on 
the dimension d of the system, ^{d = 0) = 0, /i = 1 at the 
upper critical dimension. The dashed line is obtained from 
Eq. (|l^), i. e. by the expansion from the lower critical di- 
mension, the dash-dotted line is the expansion [|l^ from the 
upper critical dimension. The points 1, 2, and 3 correspond 
to the Bak-Sneppen model at ID, 2D, and 3D. The point la 
corresponds to the ID anisotropic Bak-Sneppen model. 
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